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Abstract 

We study an inflationary scenario with a vector impurity. We show that the universe undergoes 
anisotropic inflationary expansion due to a preferred direction determined by the vector. Using the 
slow-roll approximation, we find a formula to determine anisotropy of the inflationary universe. We 
discuss possible observable predictions of this scenario. In particular, it is stressed that primordial 
gravitational waves can be induced from curvature perturbations. Hence, even in low scale inflation, 
a sizable amount of primordial gravitational waves may be produced during inflation. 
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I. INTRODUCTION 



It is often mentioned that cosmology has entered into a new stage, so-called precision 
cosmology. Of course, it is referring to developments of observational side. From theoretical 
point of view, however, we have not yet exhausted possible phenomenology on the order 
of a few percent. Clearly, it is important to explore qualitatively new scenarios at the 
percent level. Here, above all, we would like to point out that an inflationary model with 
a few percent of anisotropy yields significant consequences. Indeed, a few percent does not 
mean the consequent effects are negligible. Rather, it provides the leading component of 
the primordial gravitational waves in low scale inflationary models which are preferred by 
recent model construction in string theory 

One may feel that to seek the anisotropic inflation is against for the basics of the infla- 
tionary scenario. Actually, the isotropy is the most robust property of inflation because of 
the cosmic no-hair theorem on the isotropization of Bianchi universes [2]. However, it is 
possible to evade the cosmic no-hair theorem by incorporating the Kalb-Ramond action [s| 
or considering higher curvature theories of gravity % |5j. In spite of the possibility, no 
one has attempted to construct any inflation models based on these ideas. The apparent 
other possibility is to break the Lorentz invariance by introducing a condensation of vec- 

n q 

tor field. The model proposed in |6J, |7| seems to be successful, however, it is known to 
be metastable , llOj. A more natural possibility would be to realize a slow- roll phase 
of vector fields like as inflaton fields in chaotic inflationary scenarios. So far, it has been 
believed that it is difficult to make the vector field slow- roll without fine tuning ll| 36]. 
Very recently, the situation has changed b y t he discovery of slow-roll mechanism for the 



vector field due to non-minimally coupling 



131 ] . Hence, an apparent difficulty to construct 



the anisotropic inflationary scenario has been resolved. At this point, it is important to 
realize that both scalar and vector fields are ingredients of fundamental particle physics 
models. Therefore, it is natural to consider both scalar and vector fields exist during in- 
flation. Of course, the vector fields should be subdominant in the inflationary dynamics 
in order to reconcile the scenario with current observational data [ijj]. In this sense, the 
vector field should be regarded as impurity. Nevertheless, the effect of the vector impurity 
on observables should not be overlooked under the current precision cosmology. 

In this paper, we propose an anisotropic inflation model with the vector impurity. To the 
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best of our knowledge, this is the first concrete model which realizes anisotropic inflation, 
exits successfully to the isotropic standard universe, and provides a framework to discuss 
interesting phenomenology. We argue that the anisotropic inflation yields the statistical 
anisotropy in fluctuations. More importantly, as is expected, the primordial gravitational 
waves could be induced from curvature perturbations through the anisotropic background. 
Hence, we can expect the correlation between the curvature perturbations and the gravita- 
tional waves. In addition to these, we point out that linear polarization of the gravitational 
waves is created, which should be observed through CMB or direct interferometer observa- 
tions. 

The organization of this paper is as follows. In section II, we present the model for the 
anisotropic inflation. In section III, we analyze the system numerically and show that the 
anisotropic inflation is realized successfully. In section IV, using the slow-roll approxima- 
tion, we obtain degrees of the anisotropy. In section V, we discuss phenomenology of the 
anisotropic inflation induced by the vector impurity. There, we emphasize that the gravita- 
tional waves can be produced through the anisotropy of the spacetime. The final section is 
devoted to the discussion. 



II. THE MODEL 



In this section, we derive the basic equations for the anisotropic inflationar y s cenario 



16( | where 



induced by the vector impurity. The set up is similar to the previous works 
scalar and vector fields are considered. However, backreaction of the vector field on the 
background dynamics is neglected there. Here, we consider the backreaction and find it 
makes significant differences. 

We consider the following action for the background gravitational field, the scalar field 
and the non- minimally coupled massive vector field A^. 



S 



d xa 



-g 



i 



— R - - (d u 6Y - V(d>) - -F^F^ - - m' 



where g is the determinant of the metric, R is the Ricci scalar, V{4>) is the scalar potential, 
m is the mass of the vector field, and we have defined F^ u = d fJj A u — d u A fJj . The equation of 



motion for A^ from above action: 



7=3, ( v^F- ) = (m 2 - f) A v 



(2) 



reduces Ao(t) = in the case of v = because of antisymmetry of F^ v . As we find the 
vector field has only spatial components, we take x-axis in the direction of the vector, 



A li = {0, A x (t), 0, 



= m . 



Now, we will take the metric to be homogeneous but anisotropic Bianchi type-I, i.e. 



rfs 2 = _^( t )2 dt 2 + e 2a(t) f e -Aa +{ t) dx 2 + ^{t) ( e ^V-{t) dy 2 + £ -2 V3a. (t) ^ 



(3) 



(4) 



where N(t) is the lapse function. With this ansatz, the background action becomes 

1 
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d 4 x-^- r e 3a 

Af 



l(-a 2 + a 2 + + al)+U 2 -Af 2 V 



+\ (X - 2a + x) 2 - ^M 2 X 2 + f \a\ + \a 2 _ - 2aa + ) X 2 



(5) 



where we have introduced a new variable X = exp(— a + 2a + )A x and defined • = dt- In 
general, this procedure leads to incorrect results. However, for homogeneous Bianchi models, 
this reduction procedure is valid as can be checked directly. 

The variational equations of motion with respect to Af, 0, X, er_, a + and a then become 
(after setting Af = 1): 



- 2 {-a 2 + al + al) + l -^ + V 



K 



^ 2 {X-2a + X) + -X 2 

3a<j) + V <j} = 0, 



— ex , h — a 
2 + 2 ~ 



2a&+ X 2 = 



X + 3oX + (m 2 - 2a + - 2a&+ - 5a 2 - a 2 _) X = 



.3a 



„3a 



^ + X 2 a- 



+ 5X 2 ) a + - 2XX - 2oX 2 



K 
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a + 3d 2 - re 2 V + -A€ 2 cr+XX + k 2 -er + + a& + m 2 ] X" 
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(6) 
(7) 
(8) 
(9) 

(10) 
(11) 



where = Notice that the effective mass squared for X, i.e. m 2 — 2a + — 2a& + — 5a 2 — a 2 
in Eq. ([8]) might look different from the mass term of X in the action ([5]). However, we can 
identify them by transforming the action ([5]) into the canonical form for X. 

Note also that if the effective mass squared of X in Eq. ([8]) is negative, the system will 
be tachyonic. To avoid this, we require the effective mass squared is positive, m 2 — 2a + — 



2a& + — 5&+ — a 2 > 0, that is, m 7^ and &± has to be sufficiently small. As we will see 
later in Eq. (]17j) . this latter condition implies the amplitude of the vector field X should 
be small. In this sense, the vector field is a kind of impurity. In the limit, kX — > 0, the 
above set of equations reduce to those of conventional inflation scenarios. Also, taking look 
at Eqs.(foT)- (llll) . we see the deviation from the conventional slow-roll dynamics comes in of 
the order X 2 . 

Although we have assumed the direction of the vector field does not change in time and 
considered only x-component of it for simplicity in Eq. ([3]), we could also incorporate other 
spatial components of the vector field. In that case, each component of the vector field 
obeys the similar equation as X. The vector field turn out to start to rotate under arbitrary 
initial conditions. If the slow-roll conditions are satisfied for each, the vector will spiral down 
slowly. It is interesting to study the consequence of this complicated dynamics. However, 
we leave it for future work. 



III. ANISOTROPIC INFLATION 



In this section, we consider the evolution of the background spacetime driven by the 
scalar field 0, which we refer to as inflaton below, in the presence of the vector impurity. 

Let us consider the universe after a sufficient expansion, a — > 00. It is straightforward in 
this limit to integrate Eqs. (Q and ( flOl) to find 



*- = °< ff + = 67^T^( 2 - YX ' + 2aA ' 2 )' (12) 

We find that the anisotropy in y-z plane, &-, will disappear. This is reminiscent of the 
comic no-hair theorem by Wald [2j. However, as long as the vector field exists 1^0, the 
anisotropy in x direction, & + , will still remain even if the universe undergoes a period of 
inflation. In the proof of the cosmic no-hair theorem for Bianchi models, the strong and 
dominant energy conditions are assumed. The reason why the anisotropy does not disappear 
in our model is that the non-minimal coupling breaks the strong energy condition. Thus, 
we can restrict our metric to the following form: 

ds > = _ df 2 + e 2a(t) [ e -^ + (t) dx 2 + e 2* + (t) ^ y 2 + ^ ] ^ 

The above result might be generalized to the modified version of the cosmic no-hair theorem. 
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FIG. 1: The phase flow in X-X plane is depicted. For various initial conditions with small ampli- 
tude of X, we have plotted the trajectories. Every trajectory converges to the slow roll attractor. 

It is interesting to examine other Bianchi type models in the presence of the vector impurity. 
In particular, it is intriguing to prove the cosmic no-hair theorem for these cases. 

Now, we numerically solve Eqs. by setting er_ = 0. We will use new variables 

H = at and X = & + below. We take V = l/2/i 2 2 as the potential for the infiaton <fi. The 
parameter of the system is the ratio m/fi. For this calculation, we set k = 1, \i — 10~ 5 and 
m = 2y / 2 x 1CT 5 . And, we took the initial values 0o = 10 and H = 3 x 10~ 4 for all figures. 

In Fig. [H we depicted the phase flow in X-X plane. We set the initial value So = and 
determined <fio using the constraint equation We see that slow- roll phase of the vector 
field is an attractor when the amplitude of the vector field is sufficiently small. For the 
appropriate parameter m/fi, the small value such as X 2 ~ 0.1 is easily attainable. Both the 
mass and the non-minimal coupling play important roles for realizing the slow-roll. They 
affect the coefficient of X in Eq. (8). In detail, the mass is important for the system to be 
stable and the non-minimal coupling cancels out the term a 2 in the effective mass squared 
of X which enable the vector field to slow-roll. Thus, both mass and non-minimal coupling 
help the vector field to slow-roll. Thanks to the slow-roll for the vector field, we have enough 
anisotropic era. After the vector field goes to the minimum of its potential, the anisotropy 
disappears. The question of when the anisotropy disappears depends on the parameter mj \x 
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FIG. 2: The phase flow in cft-cft plane is depicted. For the same initial conditions as Fig.l, we have 
plotted the trajectories. Every trajectory is degenerated irrespective of the initial conditions for X. 
As we can see, the inflaton rolls down the potential slowly The inflaton begins to oscillate around 
the minimum of the potential and the inflation ends with the reheating. Here we have plotted a 
part of the graph for the sake of visualization , though we started from (f>o = 10. 
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FIG. 3: The Hubble parameter H is plotted as a function of e-folding number N . We see a slow-roll 
phase clearly. We have e-folding number N ~ 40 for this case. 
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FIG. 4: The ratio Ti/H is plotted as a function of e-folding number. In spite of the rapid expansion 
of universe, the anisotropy remains sizable for some period relevant to observations. 

ete analysis of the initial conditions is possible as 



and the initial amplitude of X. The comp 
is done in the case of pure vector models 

In Fig. [21 the phase flow in <fi - <fi is plotted with the same initial conditions as Fig. [TJ 
The trajectories are almost same irrespective of the initial conditions for X, and show the 
slow-roll phase. After slow-rolling, the inflaton field gets its damped oscillations. Therefore, 
the anisotropic inflation ends with reheating as in the standard inflation. 

In Fig. [31 we have plotted H as a function of e-folding number iV = a — «o- Actually, 
we plotted four lines for different initial conditions S /i?o = 0, 0.05, 0.1, 0.15 with fixed 
initial conditions X = 0.5 and X = 1.5 x 10~ 5 . In the Fig. [31 all lines are degenerated. 
Irrespective of the initial conditions, we see the slow-roll phase from this figure. As well as 
Fig. [21 qualitative behavior does not depend on the initial conditions for X as long as its 
amplitude is small. 

In FigJH we have plotted E/i? as a function of the e-folding number with the same 
initial conditions as Fig. [31 We see the anisotropy remains sizable for some period and then 
dumped down to zero around N ~ 20 irrespective of its initial condition. Duration of this 
phase depends on the mass of the vector field and other parameters. We find larger initial 
condition for E makes smaller anisotropy later. This can be understood from Eqs. ([8]) and 
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(|T2|) . Suppose the initial anisotropy is larger than the value calculated from Eq. (|T2|) using 
the initial data for X,X,H, then the anisotropy has to decrease rapidly to adjust to the 
value of the attractor. During this phase, the effective mass in Eq. ([S]) becomes large due 
to the term — a + . Then, the amplitude of the vector X before entering the slow- roll phase 
decreases more than the case of smaller initial condition. From Eq. (fT2"j) . the small amplitude 
of X implies the small anisotropy during the slow-roll phase. Thus, we have understood the 
behavior in Fig. HI In the case of the smaller initial anisotropy, the opposite happens. 



Previous 
siderec 
works 21 



* t 



ie conventional inflationary scenario in Bianchi type-I model has been con- 



20[. In those cases, the anisotropy decays exponentially fast. In the recent 



the anisotropic universe has been investigated in the context of dark energy. 



It would be interesting to apply our framework to the dark energy problem. 



IV. SLOW-ROLL APPROXIMATION 

Now we consider the regime where both of the vector and the scalar field are slow-rolling. 
Thanks to the slow-roll for the vector field which breaks strong energy condition, we have 
enough anisotropic era. For the reason explained in the previous section, we can exclude 
er_ from the consideration. In the slow- roll approximation, we can write Eqs. (jH])-(jBD and 
(HDD-dull as 

K 2 

3{-H 2 + £ 2 ) + k 2 V + 1k 2 Y?X 2 + — (m 2 + S 2 - AHT) X 2 = , (14) 

3H<f> + V+ = , (15) 

3HX + (m 2 - 2HY* - 5S 2 ) X = , (16) 
2X 2 

s = W^Tvo H - (17) 

K 2 

3H 2 = k 2 V + — (m 2 - 2HT) X 2 , (18) 

where we have ignored all second derivatives with respect to t and square terms of the 
first derivatives with respect to t except for H 2 , E 2 and Since we are interested in 

the situation that the anisotropy is larger than the slow-roll parameter, we have kept the 
anisotropy S 2 in the above equations. Adding Eqs. (I14p and ([TBI , we obtain Eq. (iT7|) . Hence, 
Eq. (TT4|) is redundant. 

From Eq. (|T7|) . we see X should be small in order for the anisotropy to be small as is 
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mentioned at the end of the section II. Eliminating E from Eqs. (fTTj) and (fl8|) . we can deduce 
the Friedmann equation for this model 

H 2 _ + 5X 2 / m 2 \ 

Note that if there is no scalar field, i.e. no potential, this equation tells us that H < m, no 
matter what value the vector field takes. This contradicts the condition for the slow-roll: 
H 3> m. Thus we find that the inflation does not occur only with the vector field. 
From Eq. 071) and (ITS]) , we find 



VA. ,20, 



7(6/^ + 5X 2 )(6//t 2 + X 2 )(3/«; 2 + 2X 2 ) V 2 

Here we have taken the positive sign because H > 0. Inserting these results into Eqs. (TT5"j) 
and (fTEj) . we obtain a set of differential equations which determines and X. 

Before proceeding to next section, we would like to clarify the scenario we envisage here. 
We are supposing the existence of inflaton and vector impurity. The inflation is driven by 
the inflaton. In the case of chaotic inflation, initial conditions for the inflaton are arbitrary 
as long as the energy of the inflaton does not exceed the Planck energy. The larger initial 
amplitude for the inflaton gives the longer inflationary period. One of the initial conditions 
should correspond to our present universe. Similarly, initial conditions for the vector field 
are arbitrary. However, as we have seen in section II, the large amplitude of X induces the 
tachyonic instability. Hence, this does not correspond to our universe. The criterion for the 
stability condition 2HE <C m 2 in the slow-roll phase gives 

m 3 

• 21 

\i kX 

where we used E/H ~ k 2 X 2 /3 and H 2 ~ /t 2 /i 2 2 /6 from Eqs. (ITT]) and (ITS]) . We find 
that the inflaton field has the upper limit but still has enough parameter range for the 
anisotropic inflationary universe. For example, the parameters such as K<po = 10, kXq = 0.5 
and m/fi ~ 2y / 2, which we used in the numerical calculations, satisfy this condition. 



V. IMPLICATIONS FOR PRIMORDIAL FLUCTUATIONS 



Although mixing between tensor and scalar per turbations in anisotropic spacetime has 
been studied in the literature (see, for example, (l9|), there has been no concrete realization 
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in the context of inflationary scenario. In this section, we would like to refresh various effects 
with our concrete anisotropic inflation model. Here, we do not show the explicit numbers. 
However, from the structure of the background equations, we guess possible important 
consequences of vector impurity. Especially, we stress the gravitational wave production 
due to the mixing in the context of low scale inflation. Detailed calculations will appear in 
the follow-up paper [23|. 

Since the expansion is anisotropic, we expect statistically anisotropic fluctuations. In 
fact, in the slow-roll phase (H, £ ~ const.) , we have the metric 

ds * = _ dt 2 + e 2Ht [ e ~^t dx 2 + e 2 Si ( dy 2 + ^2)] ^ 

Now, let us consider a test scalar field ip in this background metric. According to Ref. {(|, 
deviations from isotropy come in to the power spectrum of the form: 

P^k) = P (k) [l + 0(K 2 X 2 )(k • n) 2 ] , (23) 

where P^ is the power spectrum of the scalar field ijj, Po(k) is its isotropic part, n is the 
unit vector in the direction of x and k is the unit vector along the direction of wave number 
vector k. Here, in evaluating the deviations, we used the approximate relation E/if ~ k 2 X 2 
derived from Eq. ( 1TTI) . Hence, if Po(k) has a flat spectrum, the total spectrum should also be 
flat, even though it is anisotropic. In the case of curvaton scenario, the scalar field ip decays 
after reheating and converts to curvature perturbations. For this case, the fluctuation of 
i/j gives CMB fluctuations directly. Even for more complicated case of ours where inflaton 
also fluctuates, we expect the similar result as Eq. (1231) . This spectrum yields the anomaly 
in CMB observations. In paper [jjj], it is pointed out that about 10% deviation from the 
isotropic statistic can be detected by WMAP and 2% by PLANCK. Hence, the interesting 
number is around 0(k 2 X 2 ) ~ 0.1. 

In addition to this apparent effect, we can expect much more interesting ones. First of all, 
we should recall that the primordial gravitational waves are created quantum mechanically in 
the standard inflationary scenario. Hence, their amplitude should be of the order of H/M p i 
(M p i is the Planck mass). On the other hand, here, we have another new component of 
tensor fluctuations induced by curvature perturbations through the anisotropic background 
expansion. This can be understood from the background equation (TIT]) : 

d \ 1 2 Y 2 



, (24) 
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where we have used the fact k 2 X 2 -C 1 and assumed X is almost constant. If we define the 
following quantity which is supposed to be the gravitational wave: 

H+ = 5o + - -K 2 X 2 5a, (25) 
3 

we find this is a gauge invariant quantity in the long wavelength limit. The background 
equation fj24jl implies H + is a conserved quantity, namely, dH + /dt = 0. In the conventional 
isotropic case, this equation reduces to d5a + /dt = 0, which is equivalent to the conservation 
law for gravitational waves h + = in the long wavelength limit. Here, the polarized com- 
ponent of the transverse traceless tensor, h + , goes to Sa + in the long wavelength limit, that 
is, h + ^— ^ 5a + . Now, we have an extra term 5 a in the gauge invariant gravitational waves, 
so 5a + itself does not conserve in the long wavelength limit. Since the curvature perturba- 
tion is dominant during the inflation, the gauge invariant gravitational waves are generated 
first through 5a which is supposed to be the curvature perturbation. Thus, we have the 
relation H + = — ^K 2 X 2 5a. As the vector field goes to zero, 5a + starts to appear because of 
the conservation law. Finally 5a + becomes dominant after X vanishes. Thus, we have the 
relation H + = 5a + . Because of the conservation law, we have 5a + \f = — ^K 2 X 2 5a\k, where 
the left hand side is evaluated at the final time and the right hand side is evaluated at the 
horizon crossing time for the mode with wave number k. Now, let us define the curvature 
perturbations, TZ, which goes to Sa in the long wavelength limit, that is, TZ — > 5a. Then, 
we obtain the final result 

h+ = 0{k 2 X 2 )1Z , (26) 

where the right hand side of Eq. (|26|) are evaluated at the horizon crossing time. For the other 
polarization mode, we have the standard conservation law h~ = in the long wavelength 
limit. Here, we have the correspondence h~ <5cr_. 

The above discussion tells us the scalar perturbation TZ, which comes from the fluctuations 
of the trace part of metric ffl3l) . induces the tensor perturbation h + , which is the one of the 
traceless part of the metric. Therefore, the tensor power spectrum should be related to 
the curvature power spectrum Pr as 

P h+ = 0(n 4 X 4 )P n (27) 

where the relation (126!) is used. The other power spectrum P^- coming from <5<t_ has no cor- 
relation with the curvature perturbations because there is no conversion from the curvature 
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perturbations in this case. In other words, we expect linear polarization in the primordial 
gravitational waves. This could be detected through CMB or directly by DECIGO [3o| . 



In papers 



24 



25 



26, 



27 



28 



29j , the gravitational waves generated from the second order 
perturbations have been studied. It should be stressed that the mechanism to produce the 
gravitational waves we are discussing is the first order effect in the anisotropic inflationary 
universe. The primordial fluctuations produced by the anisotropy are not the conventional 
ones created from quantum fluctuations directly. They are induced from the curvature 
perturbations created from quantum fluctuations. Since the scalar perturbations are always 
larger than the tensor perturbations, this anisotropy induced mechanism is very efficient. 
This is an extremely important result. Because this implies the existence of the primordial 
gravitational waves even in low scale inflationary scenarios. In fact, supposing that we 
detect the anisotropy n 2 X 2 ~ 0.1, we would detect the primordial gravitational waves with 
the tensor-scalar ratio r ~ 0.01. 

Since at least part of the tensor perturbations are induced from the curvature perturba- 
tions, there exists the correlation between the curvature and the tensor perturbations 

(h+TZ) ~ 0(k 2 X 2 ) {Tin) , (28) 

where we used the relation (1261) . The correlation should give non-zero TB correlation at the 
10% level in the case of k 2 X 2 ~ 0.1. In other words, the normalized correlation function 



(TB) 



0(1) does not vanish. Surely, this should be searched observationally. 



y/(TT)(BB) 

Although it is necessary to calculate precise numbers for prediction, we have uncovered 
the new possibility to generate the primordial gravitational waves on dimensional grounds 
based on the concrete anisotropic inflationary model. All of the above set of predictions are 
peculiar to the anisotropic inflation and hence can be regarded as signatures of the vector 
impurity. 



VI. CONCLUSION 



We have explored the cosmological implications of the vector impurity. It turned out 
that the vector impurity affects the cosmic no-hair theorem. Consequently, the anisotropic 
accelerating universe is possible in the presence of the vector impurity. To prove it, we have 
numerically solved the equations and shown the phase flow from which we see the anisotropic 



13 



inflation successfully occurs and ends with oscillation. We found the explicit formula to 
determine the anisotropy of the inflationary universe by using the slow-roll approximation. 

We have also discussed possible consequences of the anisotropic inflation on the cosmo- 
logical fluctuations. Since rotational invariance is violated, the statistical isotropy of CMB 
temperature fluctuations can not be expected. It is intriguing to seek for the relation to the 



large scale anomaly discovered in CMB by WMAP [3l|, |32, |33| . More interestingly, the ten- 
sor perturbations could be induced from the curvature perturbations through the anisotropy 
of the background spacetime. One immediate consequence is the correlation between the 
curvature perturbations and the tensor perturbations. This correlation should be detected 
through the analysis of temperature-B-mode correlation in CMB. This new possibility im- 
plies, even in the low scale inflation, we can expect the primordial gravitational waves. This 
is an important result for future observational planning, because there has been worry that 
string cosmology tend to suggest the low scale inflation. Moreover, because of the anisotropy, 
there might be linear polarization in the primordial gravitational waves. This polarization 
can be detected either through the CMB observations or direct interferometer observations. 
These predictions can be checked by future observations. 
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b make these predictions more precise, we need to develop the perturbative analysis 
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34j . These are now under investigation 23[. The calculation of the perturbations 



is much more complicated due to the violation of rotational invariance. However, since the 
anisotropic inflationary universe is smoothly connected to the isotropic radiation dominant 
phase, the interpretation of the results should be clear. The implications in primordial 
magnetic fields and the structure formation of the universe [35] should be also studied in 
future work. 
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